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Abstract 



We calculate analytically the exact dynamical response of a droplet of 
interacting electrons in a quantum dot with an arbitrarily time-dependent 
parabolic confinement potential uj{t) and a perpendicular magnetic field. We 
find that, for certain frequency ranges, a sinusoidal perturbation acts like an 
attractive eff^ective interaction between electrons. In the absence of a time- 
averaged confinement potential, the electrons can bind together to form a 
stable, free-standing droplet. 

PACS numbers: 73.61.-r, 03.65. Ge, 73.23.-b 
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Quantum dots have attracted much interest recently, both from a pure and apphed 
viewpoint. According to the dot fabrication, the confinement length scales in the three 
spatial directions can be quite different yielding quasi one-, two- or three-dimensional dots 
||1|] . In addition the number of electrons in the dot can be reduced down to the few-electron 
limit. Most theoretical quantum dot research has been concerned with the linear response of 
the resulting A^-electron, low- dimensional system. External, time-dependent electric fields 
are typically treated as small perturbations which merely give rise to transitions between 
the eigenstates of the unperturbed dot. Given the technological possibilities for preparing 
stronger fields, it is interesting to consider the effects of larger, time- dependent perturbations 
on such dots. For example, a confinement potential with sinusoidal time-dependence could 
be created by applying an a.c. bias to the electrodes defining the dot in a heterostructure 
sample. Unfortunately, the quantum-mechanical problem of interacting electrons in a 
dot subject to an arbitrarily strong, time-varying perturbation is too complicated to solve 
in general, even numerically. 

Here we provide an analytically solvable model for the dynamical response of an A^ 
electron dot with an arbitrarily strong, time-varying confinement potential, in the presence 
of an arbitrarily strong magnetic field B. The analytic tractability of the model is made 
possible through a combination of a parabolic form for the dot confinement potential and 
an inverse- square electron-electron interaction potential [l/r^ with n = 2). The parabolic 
confinement is a reasonable approximation for many semiconductor quantum dot samples 
The true repulsive interaction between electrons in the dot is likely to be better fit by 
n ~ 1 at small r and n ~ 3 at large r due to image charge effects in neighboring electrodes 
0; however the general features of our results with n = 2 for all r should be qualitatively 
correct. Here we focus on the usual quasi- two-dimensional dot; however, most of our formal 
results can be generalized to both quasi-one- and quasi-three-dimensional dots at i? = 0. 

The time-dependent Schrodinger equation for the dot with a time- dependent confinement 
in the xy-plane, subject to a constant B field applied along z, is given within the effective- 
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mass approximation by iJ(t)\l/(t) = ih^'^it) with 




, position and z-component of angular momentum of 



given by pj, rj and /j. The cyclotron frequency is ujc (N.B. we can only solve with 7^ for 
a two-dimensional dot). The spin part of the Hamiltonian is time- independent and therefore 
decouples for all t. In order to include both transient and steady-state responses, we consider 
the time-dependent dot potential to be turned on at time t = 0: 



where uoq is the characteristic frequency of the parabolic confinement potential for t < 0, 
and f{t) has arbitrary functional form and magnitude. 



the essential results. Standard Jacobi coordinates are employed, Xj (i = 0, 1, . . . , — 1) 
where Xq = ;^ Z]j (center-of-mass), Xi = ^J^{r2 — ri), X2 = \/f(^^^^^^ ~i"3) etc. together 
with their conjugate momenta Pj. The center-of-mass motion decouples, H = HcmO^o) + 
iJrei({Xj>o}), hence E = Ecu + -^rei and \1/ = '?/'cM'?/'rei- The exact eigenstates ipcM of Hqm 
and eigenenergies -Ecm are identical to those of a single particle in a parabolic potential. 
The non-trivial problem is to solve the relative motion equation ifreiV^rei = -^'rei'^rei- We 
transform the relative coordinates {Xj>o} to standard hyperspherical coordinates: Xj = 
rdljL'i^ cosaj+i)sinaie*^* with r > and < < | (ai = |). Physically, the hyperradius 
r is just the root-mean-square electron-electron separation. The exact eigenstates of ifrei 
have the form ■j/'rei = R{r)E{Q) where Cl denotes the (2A^ — 3) hyperangular {6, a} variables; 
R{r) and E{Q) are solutions of the hyperradial and hyperangular equations respectively. 

For t > 0, the separation of the relative motion is still exact, i.e. = 4'cM(t)i^reiit) . 
The hyperangular equation is independent of uj{t), hence E{Q) remains time- independent. 
The time-dependence of \& is only contained in the center-of-mass and hyperradial parts. 
Consider explicitly i'reiit) and hence the hyperradial part R{r,t); the solution for i^cMit) is 
exactly analogous. Following Ref. p we construct the generating function 





for t > 



(2) 



For t < 0, co'^(t) is time-independent. The problem is treated in Ref. 0]; here we review 
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g{z,r,t) = J2R„.{r,t)z^ (3) 

n=0 

where Rn{r,t) are the solutions of the time-dependent hyperradial equation. Because the 
dot potential is constant for t < 0, we can employ the static Rn{r) from Ref. to obtain 

OO 2 _ r2 

g{z, r, t < 0) = E (^)"^r^-^(7j)e (4) 

n=0 '0 to 

where L denotes the Laguerre polynomial and /q = h{'m*uJo{B))^'^ . The parameter 7 is 
determined by the ujq, B and t -independent hyperangular equation, which does not admit 
complete exact solution. Fortunately, we do not need any further knowledge of the properties 
of 7 for the analysis in this paper. Equation (4) can be written in closed form using a known 
identity |] 

g{z,r,t< 0) = (^)^e^^^^'(l - ^)-("+^) (5) 
where a = 7 + — 2. We make the ansatz 

^?(^,r,t>0)=«(^,t)(fre"'(^'*)^' (6) 

which can be shown to satisfy the time- dependent hyperradial equation 

, 2N-3d 7(7 + 2iV-4) ,m\. , , , 2im* d , , 

+ TT - 2 - {^f^\ty)9{^. r, t) = T-j-giz, r, t) (7) 

or jjj. ^ Ti at 

and the t = boundary condition in Eq. (5), provided 

a{z,t) = [r/(t)]-(°+^)exp[2i^(t)(a + 1)](1 - z exp[2i9{t)])-^''+^'^ (8) 

and 

= - mm - z expmm-') (9) 

where T]{t) = |?7(t)|e*^'^*-' solves the c/assica/ one-dimensional oscillator 

r^lt) + fitUt) = (10) 
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with boundary conditions //(O) = 1 and r](0) = —iojQ^B). We may then expand g{z,r,t) 
using the relevant identity 0], and compare coefficients of with the defining Eq. (3) to 
obtain the desired time-dependent wavefunctions (unnormahzed) : 

i?„(r, t > 0) = \m'"''y'expmt){2n + a + 1) + ^^|l^Wn]e"^''^n(2/') (H) 

where y = Solving Eq. (10) hence provides a complete description for the evolution 

of the initially stationary hyperradial state i?„ for t > 0. Equation (11) together with the 

exactly analogous expression for the center-of-mass entirely determines the time-evolution 

of the total wavefunction 

First we consider a step-like perturbation (see Fig. 1(a)) 

(u^ for t < 

to^t) = \ 2 2 . (12) 

\m=^l-^ fort>0 

Increasing Q reduces the time period over which the change occurs. The solution to Eq. 

(10) for t > is 

r]{t) = {nt + l)-^{AJ,[tu^{t + fi-i)] + BN,[tu^{t + n-^)]} (13) 

where Jy and are Bessel functions of the first and second kind respectively, v = ( "^2 ° + 
^)5, and A and B are complex constants chosen to satisfy r](0) = 1 and 7^(0) = —itOQ. This 
complex function is displayed in Fig. 1(b), the inset on the right showing the view when 
projected on the complex 77 plane. For t < the value of r]{t) lies on the unit circle in the 
complex plane, and describes a simple helix around the time axis as t varies. For t > 
the function describes an ellipsoidal path in the complex 77 plane, the ellipsoid lying entirely 
within the unit circle; the eccentricity of the ellipsoid is determined by the sharpness of 
the transition from ljq to Uoo, i-e. by Q. Since Rn{r,t) (Eq. 11) depends on y = -j^, 
the ellipsoidal path of //(t) implies that the total wavefunction will oscillate for all t > 0. 
Figure 1(c) shows the time-dependence of the expectation value of the hyperradius f (r.m.s. 
electron-electron separation). Remarkably, this quantity oscillates for all t > regardless of 
the values of ujq, uJoo and Vt. Changing Vt alters the oscillation amplitude but not the period. 



In order to obtain a large amplitude of oscillation, must be of comparable order to ujq'^; 

if hujQ ~ 1 meV, should be of the order of picoseconds. 

Second we consider a sinusoidal perturbation (see Fig. 2(a)): 

( uj^ for t < 

u;\t) = \ (14) 
[ f{t) =ujI- ujI{1 - cos{2Qt)) for t > 

The solution to Eq. (10) for t > is a Mathieu function: 

oo oo 

r]{t) = Ae^™ ^ C2„e2^"™ + Ee"'^™ ^ Ca.e-^^""* (15) 

n=— oo n=— oo 

where /i and {c2„} are determined by simultaneous equations (see in Ref. 0]), and A and 
B are complex constants chosen to satisfy 7^(0) = 1 and //(O) = —iujo. Figure 3 shows fi as 
a function of ujq, uji and Q. The circle in Fig. 3 indicates the parameters for Fig. 2(a). In 
certain regions (white) of the parameter space fi is purely imaginary while in others (dark) 
it has a real part. This real part indicates that the classical particle is resonating with the 
oscillating dot; the particle's oscillations become infinitely large as t ^ oo. The correspond- 
ing effect on the quantum mechanical system, which depends on time only through ri{t), 
will be an increase in energy as t — ^ oo and a decrease in localization of R{r, t) and hence 
iprei- This spreading in -^rei implies an increase in the average electron-electron separation, 
and will lead to electrons escaping from any realistic dot having a finite depth. For small 
uji, fi has a real component only when the perturbing frequency Q is equal to an integer 
fraction ^ of the dot potential c^o- Note our treatment is exact for any amplitude uf. The 
theory remains valid even when cuf > ^Uq, in which case the dot becomes repulsive (i.e. 
an 'anti-dot') for part of each oscillation as shown in Fig 2(a). The region of Fig. 3 below 
the line u = ui corresponds to such a system. Stable zones persist here, thus a dot may 
periodically become repulsive and yet confine its electrons for all time. Indeed in the limit 
uJi = ujq (the bottom edge in Fig. 3) the N -electron droplet may remain localized for all 
time even though the time-averaged confinement potential is zero (i.e. o)^ = cug — tu^ = 0). 
The sinusoidal perturbation in these white, stable zones generates an effective attractive 
interaction between electrons which competes with the electrostatic repulsion. At = the 

6 



result is a stable, free-standing A^-electron plasma droplet. Prom the inset in Fig. 3 we see 
that, for a; = 0, the range of the parameter ^ over which the system is stable (i.e. white 
zones) is small but finite for larger we note that there an infinite number of such stable 
zones which are distributed along the entire u) — axis. 

Varying only Q corresponds to moving along a radial line in Fig. 3. Consider the 
indicated line uj = cui, moving out along this line switches the system between unstable 
states (dark zones - system absorbs energy from the perturbing field) and stable states (white 
zones - no net absorption). This pattern seems to persist arbitrarily far along this radial 
line, hence an experiment to study this stability effect could be performed at any frequency 
Q which is convenient, irrespective of cuq. This observation is encouraging experimentally 
since the frequency loq is typically high ( lO-*^^ Hz for tiwo — ImeV). 

We now focus on 'stable' solutions (white zones in Fig. 3). Figure 2(b) shows the time- 
variation of the r.m.s. electron-electron separation f{t) for a sinusoidal perturbation (Fig. 
2(a)); it seems to show a chaotic irregularity (N.B. the amplitude is always finite since we 
are in the stable regime). The function f{t) is actually the square-root of a sum of cosine 
functions; Fig. 2(c) shows the amplitudes and frequencies of the dominant terms in this 
sum. They are made up of two distinct groups: one at frequencies 2nft and another at 
(2n -|- l)fl ± pn. Two lines from the second group are labeled as an example. The value 
of (3 may be found from Fig. 3; for the present choice of parameters (marked by a circle 
in Fig. 3) /3 — 2.35. Figure 3 shows that when the perturbation amplitude cui is small 
then /? ~ ^. Thus for weak-to-moderatc perturbations, the second group of peaks in the 
frequency spectrum lie at sums and differences of Cu and the driving frequency il, i.e. at 
(2n -|- 1)Q ± (D. The system can hence exhibit significant frequency mixing and harmonic 
generation. Figure 3 shows how (3 (and hence the spectrum) changes as we move to the 
strongly non-linear regime {cui ~ Cu). The parameters for Fig. 2 correspond to a strong 
perturbation {cof — 0.57a;o) which causes fip to deviate from O by 16%. 

In summary, we have presented a method for finding the exact dynamical response of an 
electron droplet. Our first example shows that a droplet subjected to a sudden increase in 
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confinement will exhibit novel size oscillations associated with dynamical breathing modes 
for all time. Our second example shows that a sinusoidally varying confinement generates 
strongly non-linear size oscillations whose stability depends non-trivially on the confinement 
parameters. Electrons may remain localized even in the absence of a time-average confine- 
ment potential, thereby offering a possible new method of quantum dot fabrication. This 
plasma droplet could be moved around in the plane under the infiuence of a small, static 
in-plane electric field. In both examples the size oscillations are strong in that they involve 
all the electrons in the dot, and should therefore be readily observable experimentally. 

We thank Luis Quiroga and Nikos Nicopoulos for useful discussions. S.C.B. is supported 
by an EPSRC studentship. 
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Figure Captions 

Figure 1. (a) Step-like perturbation u;'^{t). (b) r]{t); inset shows the projection onto the 
complex T) plane, (c) The r.m.s. electron-electron separation f(t). 

Figure 2. (a): Sinusoidal perturbation u'^{t); = Uq —uj. (b) The r.m.s. electron-electron 
separation f(t). (c) Frequency (Fourier) composition of f^(i). 

Figure 3. Top: Contour plot showing as a function of perturbation parameters; /j, — 
i?e{/x} -I- ip. White regions: i?e{//} = and contours are lines of constant /3; on a;i = axis 
P — fi- Dark regions: Re{n} > and contours are lines of constant i?e{/x} in increments of 
0.2, as shown by key. 
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